function result = square(line,sav_path,info)
    global T;
    len = length(line)^(1/2); 
    result = zeros(len);
    for i=1:len
        l = len + 1 - i;
        ll = line((i-1)*len+1:1:i*len);
        result(l,:) = ll;
    end
    
    fid = fopen(sav_path,'a');
    fprintf(fid,[info,'\n']);
    for i=1:len
        fprintf(fid,'%7d',i-(T+1));
    end
    fprintf(fid,'\n');
    for i=1:len
        fprintf(fid,'%-3d', (T+1)-i);
        for j=1:len
            fprintf(fid,'%7.4f',result(i,j));
        end
        fprintf(fid,'\n');
    end
    fprintf(fid,'\n');
    fclose(fid);
end
